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Abstract 

At high level density, two states avoid usually crossing at the critical value a CT of the parameter 
a by which the system is controlled. The wavefunctions of the two states are mixed in a finite 
parameter range around a CT . This holds true for discrete states as well as for narrow resonance 
states which are coupled via the environment of scattering wavefunctions. We study the influence 
of avoided level crossings onto four overlapping complex eigenvalues of a symmetric non-Hermitian 
operator. The mixing of the two wavefunctions around a cr is simulated, in each case, by assuming 
a Gaussian distribution around a cr . At high level density, the Gaussian distributions related 
to avoided crossings of different levels may overlap. Here, new effects arise, especially from the 
imaginary part of the coupling term via the environment. The results show, moreover, the influence 
of symmetries onto the multi-level avoided crossing phenomenon. 
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I. INTRODUCTION 



Avoided level crossings play an important role in quantum mechanics. They are known 
for about 80 years [1]. In the case that the Hamilton operator of the system is Hermitian 
and the real eigenvalues provide the energies of the discrete states of the system, two of 
its states never will cross. Instead they avoid crossing when controlled by a parameter. This 
phenomenon is an observable effect: it consists in an exchange of the two states, including 
their populations, in a many-particle system. The exchange takes place at the critical value 
a cr of the parameter a. The corresponding crossing points of the energy trajectories can be 
found by analytical continuation into the continuum [2j. 

Narrow resonance states show a similar behavior. They are described well by the com- 
plex eigenvalues Ei — i/2Ti and biorthogonal eigenfunctions $j of a non- Hermitian Hamilton 
operator [3]. The eigenvalues provide not only the energies E{ but also the widths I\ (in- 
verse lifetimes) of the states. Mostly, the resonance states avoid crossing similar to discrete 
states. However, in contrast to the eigenvalue trajectories of discrete states, the eigenvalue 
trajectories of resonance states can cross. The crossing points are called usually exceptional 
points according to Kato |4j who studied first the nontrivial mathematical properties of these 
singular points. Also in the case of avoided crossings of resonance states, the corresponding 
crossing points can be found by varying one additional parameter. 

Recent studies have shown qualitatively that dynamical phase transitions in open quan- 
tum systems occur when the level density is high and, correspondingly, many neighbored 
states avoid crossing in a relatively small parameter range As a result, the system is 
(dynamically) stabilized: the number of states is reduced (by one when the environment con- 
sists of one continuum), the individual spectroscopic properties of the original states are lost, 
and the narrow (trapped) resonance states of the system show collective features. Mean- 
while, unexpected experimental results from different fields of physics could be explained 
qualitatively by means of this phenomenon (see where some of them are sketched). The 
dynamical stabilization of the system is environmentally induced and appears only in open 
quantum systems which are described by a non-Hermitian Hamilton operator due to the 
embedding of the system into the environment of scattering wavef unctions. 

It is the aim of the present paper to study in detail generic features of open quantum 
systems at high level density where more than two states avoid crossing and the ratio of 
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the widths of the states to the energy differences between them is larger than 1. A realistic 



case of such a type is considered first in nuclear physics 



6] , then in laser induced continuum 



structures in atoms [7], and later in many other systems, see the review [sj and also j^]. 
We restrict the study to four overlapping resonance states and three different avoided level 
crossings. In difference to the paper jg] on the signatures of symmetric Hamiltonians with 
three coalescing eigenfunctions, we are interested in the behavior of the eigenvalues of a 
symmetric non-Hermitian operator in the regime of overlapping resonances where the states 
avoid crossing. Here, the wavefunctions of the two states are mixed in a finite parameter 
range around the critical value a cr of the avoided crossing. The mixing range shrinks to 
one point when analytically continued up to the exceptional point by means of another 
parameter {9]. At the exceptional point, the two eigenfunctions become linearly dependent 
from one another jjj. 

Our calculations are performed with real, complex as well as with imaginary coupling 

coefficients of the states via the environment. The mixing of the wavefunctions caused by 

Jl 

this coupling is simulated, according to the numerical results obtained in [9j, by assuming a 
Gaussian distribution around the critical point of avoided crossing at the critical parameter 
value a cr . At high level density, the Gaussian distributions related to avoided crossings of 
different levels may overlap. Our results show that, in the regime of overlapping Gaussian 
distributions, new effects arise from the imaginary part of the coupling term via the en- 
vironment which are related to width bifurcation. They are much more pronounced than 
the effects caused by the real part of the coupling term which leads to level repulsion in 
energy. Furthermore, symmetries between the states are shown to play an important role 
in the avoided level crossing phenomenon. This result hints at the relation of avoided level 
crossings to exceptional points. 

In section [Til t ne formalism used for the calculations, is formulated. Some results for two 
crossing levels are given in section II III while those for four crossing levels can be found in 
section IIV1 The results are summarized and conclusions are drawn in the last section. 
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II. FORMALISM 



We consider an N x N matrix 



U 



^21 ^2 



(1) 



y WATl U}]\J2 ■ ■ ■ £n J 

the diagonal elements of which are the N complex eigenvalues S{ = — z/2 7$ of a non- 
Hermitian operator. The and 7, denote the energies and widths, respectively, of the N 
states without account of the interaction of the different states via an environment. This 
interaction is contained in the values which stand for the coupling matrix elements 
{(j)i\H\(j)k) of the states % and k via the environment, where <pi is the wavefunction of the 
state i. The corrections due to the coupling uja of the state i to the environment (i.e. to 
the continuum of scattering wavefunctions into which the system is embedded) lead to the 
selfenergy of the state i, see jsj. In atomic physics, these corrections are known as Lamb 
shift. With the only exception of figure [3 the Ua are assumed, in our model calculations, 
to be included into the £j. The ui^ are complex, generally. The border case of purely 
imaginary corresponds to frozen internal degrees of freedom For discrete states with 
real energies Ei = e^, and the Uik are real {3, 

The eigenvalues of % will be denoted by Si = E{ — i/2 T{ where Ei and Tj stand for 
the energy and width, respectively, of the eigenstate i. The eigenfunctions $j of H can be 
represented in the set of basic wavefunctions <pi of the unperturbed matrix (corresponding 
to the case with vanishing coupling matrix elements o;^-), 



N 



(2) 



The bij are normalized according to the biorthogonality relations of the wavefunctions 



(3) 



It follows 



($,1$,) = Re (($,!$,)); ^ = <$ f |$ 4 )>l 



(4) 
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and 



($^> = z Im = -<*i*|*i> ! 1^1 = m\^\ > . (5) 

The £i and $j contain global features that are caused by many-body forces induced by the 
coupling of the states via the environment, see 0, Q] and equations (7) and (8) in 

Some years ago, the case N = 2 with = ej(a), fixed real w = u>u = 0J21, and different 
fixed values of including 7, = 0, is studied as a function of the parameter a in the 
neighborhood of avoided and true crossings of the two levels In this case, the two 
eigenvalues of H are 



e i,j = e i,j - = £l+ 2 £2 ±Z '' Z = \^ £l ~ £2 ^ + 4uj2 ■ ^ 

According to this expression, two interacting discrete states (with 7^ = 0) avoid always 
crossing since u and E\ — 62 are real in this case 3[. Resonance states with nonvanishing 
widths 7j repel each other in energy according to Re(Z) > while the widths bifurcate 
according to lm(Z). The two states cross when Z = 0. 

CI 

The results for the N = 2 case [9| show further that the wavefunctions of the two states 
$1 and $2 are mixed in a finite range of the parameter a around the critical value a cr at 
which the two states avoid crossing. This holds true not only for resonance states but also for 
discrete states. Furthermore, a nonlinear source term appears in the Schrodinger equation 
in the neighborhood of an exceptional point and the critical point of an avoided crossing, 
respectively. This source term causes irreversible processes 

In our calculations, the mixing coefficients bij (see equation (J2J)) of the wavefunctions of 
the two states due to their avoided crossing are not calculated. We simulate the fact that 
the two wavefunctions are mixed in a finite parameter range around the critical value of 



their avoid crossing [9] by assuming a Gaussian distribution 

co^j = lo e~^-^ 2 (7) 

for the coupling coefficients. 

Calculations for iV > 2 and for the case with complex u are not performed up to now. Of 
special interest is the situation at high level density where the ranges of avoided crossings, 
defined by jTJ, of different levels overlap. Some generic results obtained with 2 and 4 
resonance states will be presented in the following sections. 



III. N = 2 CROSSING LEVELS 




FIG. 1: The energies (top) and widths Tj/2 (bottom) for two crossing states as a function of 
the parameter a with e\ = 1 — a/2; = a; 7i/2 = 0.5; 72/2 = 0.5999; uj = 0.05. The dashed 
lines show the unperturbed ej(a). The exceptional point is at the crossing point of the energies e\ 
and e.i- 

First we study the two-level case most properties of which can be found in the literature 
at different places. Here, we choose the matrix ([I]) with = ej(a), different fixed values of 
7i and fixed u = U12 = 0021- The functional dependence of the energies over the parameter 
a is similar as in [9|, but the u may be complex. We show the results as a function of the 
parameter a in the neighborhood of avoided and true crossings of the two levels. 

In figure [U the energies Ei and widths are shown for two crossing levels with the real 
coupling coefficient u> = 0.05. At the critical value a cr of the parameter a, the two energy 
trajectories ej(a); i — 1,2 cross. Here, also the energy and width trajectories E^a) and 
Tj(a), respectively, may cross as shown in figure [TJ The crossing point of the Ei(a) and 
Tj(a) is an exceptional point. If one of the input values is slightly different from the values 
used in figure [U the Ei(a) avoid crossing while the Tj(a) cross (or the Ei(a) cross and the 
widths Ti(a) bifurcate). In the first case, the two states are exchanged at the critical point 
a cr of the parameter a. 

The avoided crossing phenomenon of discrete states is known in the literature for many 
years As shown in [9|, narrow resonance states avoid crossing in energy in a similar 




FIG. 2: The energies £7j (top) and widths Tj/2 (bottom) for two crossing states as a function of 
the parameter a with e± = 1 — a/2; e2 = a; 7i/2 = 0.5; 72/2 = 0.5980; oj = 0.05(1 + i). 
The dashed lines show the unperturbed ej(a). The exceptional point is shifted in a relative to the 
crossing point of the energies e\ and e^- 

manner as discrete states. The only difference to the case with discrete states are the 
nonvanishing widths Tj(a) of the resonance states. The widths may cross at the critical 
parameter value a = a CT if u is real and the E^a) avoid crossing, see equation (J6]). 

In figure El the situation with two crossing levels but complex coupling u is shown. Also 
in this case an exceptional point can be found, though for a value of 72 being different from 
that in figure [1] due to the other value of u. In figure [2j the position of the exceptional 
point is shifted in relation to the crossing point of the ej(a) trajectories. Furthermore, the 
difference |Ti — T 2 | blows up around the crossing point of the e^. However, at large distances 
of the parameter a from the critical value, the widths of the two states approach quickly 
the values 71 and 72, respectively. If one of the input values is slightly changed, the two 
eigenvalue trajectories avoid crossing in energy similar as in the case with real coupling uj 
(figured]) while the Tj(a) cross freely, or the Ei(a) cross and the Tj(a) bifurcate. In the first 
case, the two states are exchanged at the critical parameter value at which the two levels 
avoid crossing. 

The blowing up of the difference \Ti — around a CT seen in figure EJ is nothing but 
width bifurcation according to Im(Z), equation ([6]). This can be seen very clearly in figure 
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[3] calculated with imaginary u. In a relatively small parameter range around the crossing 
point of the e\ and €2 trajectories, we have E\ = E2 while the difference \Ti — ^1 is large. 
Beyond this parameter range, quickly Ei — > ej and Tj — > jj. The width bifurcation is caused 
by Im(u;). It plays an important role for resonance states that cross at high level density 
since, in difference to the avoided crossing of discrete states, the coupling of resonance states 
via the environment is complex in this case |3J]. 

Figures [1] to [3] show the characteristic features of the dynamics of an open quantum sys- 
tem in the neighborhood of an exceptional point: the two crossing states are exchanged at 
the critical value of the control parameter including their population and, furthermore, the 
widths bifurcate due to Im(o;). Population transfer occurs in realistic systems. A few ex- 



amples are the high-order harmonic generation in a driven two-level atom [TJj, the ultrafast 
stimulated Raman parallel adiabatic passage by shaped pulses [12 1 and the molecular vi- 
brational cooling by adiabatic population transfer from excited to ground vibrational states 
131 ] . Width bifurcation becomes important at high level density as will be shown in the 
next section. 
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FIG. 3: The energies Ei (top) and widths Tj/2 (bottom) for two crossing states as a function of 
the parameter a with e\ = 1 — a/2; e<i = a; 7i/2 = 72/2 = 0.5; uj = 0.05i. The dashed lines 
show the unperturbed ej(a). The widths bifurcate in the parameter range where E\ = E2. 
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FIG. 4: The energies Ei (top) and widths Tj/2 (bottom) for one state crossing three other states 
as a function of the parameter a with e\ = 1 — a/2; e 2 = 1.05 — a/2; = 1.1 — a/2; e4 = a; 
7i/2 = 0.5; 7 2 /2 = 0.4; 73/2 = 0.6; 74/2 = 0.58523; w = 0.05 (1 + i). The dashed lines show 
the unperturbed ei{a). The widths Tj of the three overlapping states 1, 3 and 4 bifurcate (or cross 
freely) in the parameter range of intersection while energy E2 ~ e 2 and width T2 ~ 72 remain 
almost unaffected. 

IV. N = 4 CROSSING LEVELS 

Some results of calculations with one resonance state crossing three other ones are shown 
in figures H] up to [91 Mostly we use a linear functional dependence of the energies ej(a) on the 
control parameter a. Figure H] shows the role, the symmetry in the overlapping of resonances 
plays in the critical region of the parameter a. Energy e 2 and width 72 of the state 2 are fully 
symmetric with correspondence to the two states 1 and 3. The avoided crossings of the state 
2 with the states 1 and 3, respectively, are disturbed therefore completely symmetrically by 
the states 3 and 1, respectively. Due to this fact, the state 2 does (almost) not take place 
in the avoided level crossing phenomenon but appears as an 'observer'. Both, the energy 
trajectories E 2 and the width trajectories T 2 follow the trajectories e 2 and 72 everywhere 
including the parameter values a around a„. Such an effect is known from a realistic case 



with 3 crossing levels in a quantum dot 14]. 



The crossings and avoided crossings in energy together with the corresponding bifurca- 
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FIG. 5: The same as ngureHbut 7l /2 = 7 2 /2 = 73/2 = 74/2 = 0.5. u = 0.05 (1 + i). All states 
participate in the interaction scenario. Far from the critical parameter range, all Tj approach the 
value 0.5. 

tions and crossings in width of the other states 1, 3 and 4 can be seen very clearly in figure 
HJ Far from the crossing region, all eigenvalue trajectories 1^; i = 1,3,4 approach those of 
the "fi^j- An exchange of the states takes place in the critical region around a cr with the 
only exception of state 2. 

The numerical symmetry of state 2 relative to the states 1 and 3 is somewhat disturbed 
when the widths 7, of all the states are equal to one another, see figure In this case, all 
four states are affected by the crossings. In the critical region, the avoided level crossings 
in energy and free crossings in width as well as the free crossings in energy and bifurcations 
in width can be seen. Far from the critical region, all eigenvalue trajectories and Tj 
approach the trajectories Cj and jj, respectively. The small shift in energy seen in figure H] 
does not appear in figure [5] due to the same values of all 7^ in the last case. This shift is an 
indicator of the residual influence of the intersection onto the eigenvalues of % far from the 
critical region with avoided level crossings. In other words, it is a signature of the different 
couplings of the states to the environment in the regime of overlapping. 

In difference to the results shown in figures H] and [51 in those of figure |6] the states 1, 2 
and 3 do not interact via the environment: OJ12 = OJ13 = U23 = 0. Their interaction with 
state 4 via the environment is = = wej(a)e _(ei_e4 ' )2 . The results are very similar to 
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FIG. 6: The same as figure[5]but oj^ = oj^ = oj e.j(a) e ~( ei_e4 ) 2 ; w$ = ojjjtA i = 0; W44 = 0. 
lo = 0.05(1 +£). The widths bifurcate less than in figure[5J All states participate in the intersection 
scenario. Far from the critical parameter range all Tj approach the value 0.5. 

those shown in figure [5j except for the fact that the overall interaction is smaller. 

Figure [7] shows the influence of the selfenergy term W44 7^ 0. It causes some shift in all 
Ei and all Tj what can be seen best far from the critical region. 

The relative influence of the imaginary part of the coupling strength u in the region 
around a CT can be seen in figure [HJ The calculations are performed with the same values as 
in figure E]but u = 0.05(1 + i/ 10) instead of u = 0.05(1 + z). The width bifurcation in figure 
IE] is smaller than that in figure |6] due to the smaller value of Im(u;). 

In figure [9j we show the results for a completely other distribution of the levels in energy. 
The functional dependence of the energies is chosen according to a Coulomb-like potential. 
All the features discussed in figures H] to [8] can be seen also with this distribution. Even the 
influence of the symmetry of the states 1 and 3 in relation to the state 2 causes the state 2 
to be an 'observer', in the same manner as discussed in figure HI Further results obtained in 
other calculations with the distribution of figure |9] are not shown in the present paper. 

Similar results are obtained also for the case that two resonance states cross two other 
ones according to, e.g., t\ = 1 — a/2; e 2 = 1.05 — a/2; e 3 = 0.05 + a; e 4 = a, and different 
values for the 7^. Numerical results have shown furthermore that the Gaussian distribution 
(I7j) used in the calculations, ensures the overlapping of the different avoided crossings, i.e. 
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FIG. 7: The same as figure[6]but CJ44 = oj = 0.05(1 + i). The widths Tj bifurcate in the parameter 
range of intersection. The selfenergy term CJ44 7^ causes shift in energy and width far from the 
critical parameter range. 
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FIG. 8: The same as figure [6] but u = 0.05 (1 + z/10). All states participate in the avoided level 
crossings. Far from the critical parameter range all Ti approach the value 0.5. 

the overlapping of the parameter ranges in which the wavefunctions of the two crossing 
states are mixed. Results with a broader Gaussian distribution are almost the same as those 
shown in the figures of the present paper. 

The features caused by the true and avoided crossings shown in the figures [H to are 
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FIG. 9: The energies E± (top) and widths 1^/2 (bottom) for four crossing states as a function of the 
parameter a with ei = 1 - ^; e 2 = 1.05 - e 3 = 1.1 - e 4 = 7i/ 2 = 0.5; 7 2 /2 = 

0.4; 73/2 = 0.6; 74/2 = 0.58523; to = 0.05 (1 +*). The dashed lines show the unperturbed ej(a). 
The widths V%, T$, bifurcate in the critical parameter range while the width T 2 ~ 72 and the 
energy E% ~ e 2 remain almost unchanged (as in figured]). 

generic. In any case, the influence of the level crossings occurring in a relatively small region 
around the critical value a CT , onto the four eigenvalues of H is, generally, large. Only at 
large distances from the critical region, the widths l -1 , of the states approach the values jj 
with j ' i in most cases. A shift occurs due to the selfenergy term. In the critical region, 
the eigenvalue trajectories are exchanged usually. Exceptions occur when the exchange is 
prevented due to certain symmetries, such that one of the states acts as an 'observer'. 

The results obtained at large distances from the critical region are true from a mathemat- 
ical point of view. They are, however, hardly realized in physical systems for the following 
reason. In the critical region, the system is split into different parts due to width bifur- 
cation: one part contains the long-lived states while another one contains the short-lived 
states. The short-lived states decay quickly, and the long-lived states form a new system the 
wavefunctions of which are strongly mixed. This process occurring in the critical region, is 
irreversible. As a consequence, the system will never reach the parameter region far from 
the critical value a cr . Instead, a dynamical phase transition takes place as discussed in 
The figures H] to [9] show that signatures of the dynamical phase transition can be seen already 
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in the case of only four overlapping resonances. 



V. CONCLUSIONS 

In the present paper, we have shown numerical results for avoided crossings of two and 
four states in a relatively large range of the control parameter a by using the matrix (pQ). 
The avoided crossing of two levels at the critical value a cr influences the complex eigenvalues 
of the non-Hermitian operator H in a small parameter range around a cr . Here, the levels 
avoid crossing in energy while the widths cross freely or they cross freely in energy while the 
widths bifurcate. In the first case, the two states are exchanged, including their populations. 
The results can be understood by means of the eigenvalue equations (J6]). 

The influence of more than one avoided level crossing occurring in a small critical param- 
eter range, onto the eigenvalues of H is stronger than in the case with only one avoided level 
crossing in the critical region. The influence can be seen still beyond the critical region. The 
imaginary part of the interaction u of the states via the environment affects the eigenval- 
ues stronger than the real part of u does. Especially width bifurcation appearing at nearby 
avoided level crossings is usually strong. The reason for these results are the nonlinear source 
terms that appear around an avoided level crossing in the Schrodinger equation {3, 5|. They 
become more important when several avoided level crossings are near to one another and the 
wavefunctions of the different states overlap. An exception from these results occurs when 
one state is related symmetrically to two neighboring states. In such a case, the state does 
(almost) not participate in the avoided level crossing phenomenon, but is solely an 'observer' 
of the avoided crossing of the two neighbored states. This result corresponds to sensivity of 
exceptional points to the distortion of their symmetry, see the discussion in appendix E in 

0. 

In physical systems u as well as 73 are usually functions of energy. Furthermore, the ratio 
Im(u;) / Re (a;) increases with energy. We did not simulate these effects in our numerical 
calculations since we are interested in the study of generic effects. In any case, the critical 
range of nearby avoided crossings affects the eigenstates of a non-Hermitian operator in a 
relatively large parameter range. Far from the critical region with overlapping resonances, 
the eigenvalue trajectories Si = — i/2 Tj approach, as it must be from a mathematical 
point of view, the trajectories Ej = ej — i/2 jj with i 7^ j as a rule. The selfenergy term 
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causes some shift of the Si relative to the £,-. 



The results shown in figures H] to |9] are relevant for physical systems at high level density 
where different avoided crossings overlap. This range is characterized by level repulsion 
and width bifurcation as the results of the present paper show. Width bifurcation limits 
the existence of the system: it causes a splitting of the system into different parts with 



significant different lifetimes. Finally, a dynamical phase transition takes place 



Further 



studies with more than four overlapping resonance states and with uj = u(a) are in progress. 



Furthermore, we will apply the results of multi-level avoidec 



atomic and semiconductor cavity QED which is studied in 
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crossing also to the case of 
Hi?]. 
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